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Abstract 

The quantity w 1 (r) = (1 / k 2 T c \)[dn(r ,T e \) / dT e \] v ^ N is introduced as a con- 
venient measure of the local isoelectronic reactivity of surfaces. It character- 
izes the local polarizability of the surface and it can be calculated easily. The 
quantity w N (r) supplements the charge transfer reactivity measured e.g. by 
the local softness to which it is closely related. We demonstrate the applicabil- 
ity and virtues of the function w N (r) for the example of hydrogen dissociation 

and adsorption on Pd (100). 
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Density-functional calculations of chemisorption processes and of potential energy sur- 
faces of the dissociation of simple molecules over surfaces, which have become available in 
the last years, have greatly improved the fundamental understanding of reactions at solid 
surfaces. However, those extensive computations are limited to a restricted number of model 
systems. It remains an important task to develop a methodology that allows the prediction 
and interpretation of reactions at surfaces in terms of the properties of the non-interacting 
systems. Such concepts, well known in molecular chemistry, are the essence of "reactivity 
theory" . They are based on low-order perturbation theory and aim at a description of the 
early stages of chemical interactions. The ensuing response functions are called "reactivity 
indices" and characterize the changes of the electronic structure of one reactant as stimu- 
lated by the presence of the other or vice versa. In an early contribution to this field Fukui 
et al. established the correlation between the frontier-orbital density, i.e., the density of 
the highest occupied and lowest unoccupied molecular orbital (HOMO and LUMO), and the 
reactivity of a system towards electron donation or accptance. Pearson H introduced the 
electronic "softness" , the magnitude of the change of the electronic structure due to a change 
of the number of electrons in the system, as a measure of the reactivity. Species are classified 
as "soft" if only a small energy is required to change their electronic configuration, i.e., if 
the valence electrons are easily distorted, polarized, removed or added. A "hard" species 
has the opposite properties holding its valence electrons more tightly @,f|. The utility of 
the hardness-softness concept is based on the so called hard and soft acid and base (HSAB) 
principle formulated by Pearson [Q which states that hard-hard (soft-soft) interactions are 
preferred. In the case of polyatomic or extended systems the HSAB principle is used in a 
local version: The soft (hard) parts of one reactant prefer to interact with soft (hard) areas 
of the other || . Parr and collaborators || gave a foundation in density functional theory to 
those mostly semi-empirical concepts, and Cohen et al. [§|-|8j have reviewed the foundations 
of reactivity theory and addressed some unresolved issues. 

The (local) softness and hardness describe the response of the electron density to a 
change of the charge state of the system. For extended, gapless systems the local softness 
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s(r) is defined as 

where v(r) is the electrostatic potential due to the nuclei, and T e \ is a temperature (see 
below). The global softness, S, is defined as the integral of the local softness over all space. 
The local softness is the local density of states at the Fermi energy, g(r,/j,), screened by a 
response function closely related to the static dielectric function 0; it describes the local 
response of the electron density to a global change of the charge state of the system. The 
chemical interaction will be, however, largely determined by the nonlocal response of the 
electron density to a local perturbation, and is governed by the whole spectrum and not just 
the states at the Fermi energy ||. Despite these restrictions, the local softness, in connection 
with the HSAB principle, has become an important predictive tool in analyzing reactions 
between molecules P/TCfl. 

^From the beginning, application of the softness concept to surface processes and, in 
particular, to catalytic reactions at metal surfaces has been proposed, and empirical rules, 
e.g. for the influence of adsorbate layers on the surface properties of metals, have been 
deduced P,[TT|-|i"3"| . However, for extended systems molecular orbitals and levels are not 
well defined and models of chemical reactions at surfaces are reformulated using a basis 



of localized atom-like orbitals and the projected density of states (see e.g. fi4|-|T7|). Only 
recently have investigations appeared which apply the HSAB principle to semiconductor 
systems |T8|-^T|. Calculations of the local softness at metal surfaces have not yet been 



done so far. Metal surfaces are characterized by a non-zero density of states (DOS) at 
the Fermi level, and screen external perturbations efficiently by low-energy electron-hole 
excitations. In the framework of reactivity theory the metal surface is thus soft. In much 
the same spirit, Feibelman and Hamann |EJ related the change of reactivity of transition 
metal surfaces in the presence of adsorbate layers to the spatial variation of the DOS at Ep. 
Although their study demonstrated that the states in the vicinity of the Fermi level govern 
the reactivity of metal surfaces, no clear correlation between local reactivity and the spatial 



variation of the DOS close to Ep could be established. Falicov and Somorjai [ 23 1 stressed the 



importance of low-energy occupation fluctuations for the large reactivity of metal surfaces. 



Yang and Parr [|HJ as well as Baekelandt et al. have shown that the local softness is 
directly related to the correlation function of changes of the total number of electrons and 
the density fluctuations. 

In this paper we introduce the change of the electron density due to the excitation of low- 
energy electron-hole pairs induced by an increased electron temperature, T e \, as a quantity 
which characterizes the spatial distribution of the reactivity of metal surfaces. We show that 
this isoelectronic reactivity index, denoted as w N (r), represents a convenient measure of the 
capability of the metal electrons to react in response to incoming atoms or molecules: 

_ n(r,T el ) - w(r,0) 

(^ei) 2 ' 1 } 

where N is the number of electrons in the system. 

The local softness and the function w (r) are closely related, as may be concluded 

already from their definitions in equ. ([!]) and equ. (|2|). Below, a formal expression for the 

function w N (r) is presented showing this connection and the physical significance of w N (r) 

more clearly. We define in addition to the isoelectronic function w N (r) a function u> M (r) by 

taking the derivative in equ. (||) at constant Fermi level \x instead of at constant number of 

electrons. Both functions are related by 

w N (r) = ^(r) -s(r)^ (4) 

where s(r) and S are the local and global softnesses and If is the integral of w^(r) over 
space. Equation |] shows that the integral of w N (r) over space is zero, as expected for an 
isoelectronic reactivity function. The density, n(r, T el ), entering Eq. 0may be expressed as 
the integral of the product of the Fermi function, /(e), and the local Kohn-Sham density 
of states g(r, e) over the energy e. In this paper we use the Fermi function at a non-zero 
auxiliary temperature, T e \. Applying the same technique as has been described for the local 



softness in Ref. 0, the function iu M (r) is the solution of an integral equation of the same 
structure as that derived for the local softness 

w »(r) = Jd 3 r' K-\r,r') i(r') . (5) 

In the limit of zero electron temperature T e \, the kernel in equ. (|5]) is the same response 
function as found for the local softness f?]. The source term i(r) 

= iJ^rfde g(v,e) (6) 

(7) 
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is proportional to the energy derivative of the local density of states at the Fermi energy. 
Thus, the function u> M (r) is essentially a (short-range) linear mapping of the energy derivative 
of the local DOS at the Fermi level. The local softness, on the other hand, is the same 
mapping but applied to the local DOS at the Fermi level J7|J|. 

As mentioned in the introduction the local softness measures the correlation between 
density fluctuations and fluctuations in the total number of electrons in a statistical en- 
semble fL2| , |24]| . In a grand-canonical ensemble the reactivity index w N (r) has an interesting 
relationship to the correlation function between the density fluctuations and the fluctuations 
of the total energy of the surface: 

™» = ^5 ^ H ~ ^» " <"(*)»> 

((H-(H))(N- (N))}} , (8) 

where H is the Hamiltonian of the system and () indicates the averaging over the grand- 
canonical ensemble. 

The quantities u> M (r) and w N (r) describe the spatial variation of the strength and the sign 
of the density change connected with low-energy excitations of the valence electrons. They 
represent essentially the difference of the densities of the LUMOs and HOMOs: Regions 
of positive values of w N (r) correlate to dominant spatial weight of unoccupied states and 



negative values to occupied states close to the Fermi level. Thus, the function w N (r) gives a 
spatially resolved picture of the essential information contained in the local density of states 
around the Fermi energy, which also forms the starting point for tight-binding-like models of 
the local reactivity |T5| - [r7J| . Moreover, the concept underlying the definition of reactivity in 



Eq. (2) is flexible; different reactivity indices can be defined incorporating weight functions 
other than the Fermi function used here for the occupation of Kohn-Sham eigenfunctions. 

The determination of w N (r) requires two self-consistent calculations of the electron den- 
sity for different electron temperatures T e \ without changing the number of electrons in the 
system, a task readily performed. In contrast, the change of the electron temperature at 
constant fi and constant external potential used defining u> M (r) results in a change of the 
number of electrons per unit cell. The resulting charged systems cannot be treated within 
conventional supercell geometries. The finite difference formulation of equ. (3) has in addi- 
tion the advantage that it enables the modelling of an increasing strength of the interaction 
of the adsorbate with the surface by increasing the value of T e \. 

As an example of the applicability of w (r) to the study of the reactivity of metal 
surfaces we consider the interaction of a hydrogen molecule with the (100) surface of Pd. 
The calculations presented below were performed using DFT together with the generalized- 
gradient approximation (GGA) [[25|] and employing the full-potential linear-augmented plane 
wave method |2^j27|1 . The parameters used in the calculation correspond to those used in 



Rcf. p8|. The metal substrate is modeled by five layers separated by a 10 A thick vacuum 



region. 

The function w N (r) is displayed in Fig. 1 for clean Pd(100), calculated using kT e \ =70 
meV and equ. (|). The redistribution of electrons is quite pronounced. There is a clearly 
identifiable electron rearrangement within the ci-shell: The (x 2 — y 2 ) electron density is 
increased (see Fig. 1 a,b) and the (3z 2 — r 2 )- (see Fig. 1 c,d) and (xy) -orbit als (see Fig. 1 
a) are depleted. The (x 2 — y 2 )-orbitals contribute to the bonding in the surface between the 
Pd atoms, which can be inferred from the fact that their band is rather wide and contains 
the highest number of unoccupied d states at the surface. 
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We now consider the dissociative adsorption of a hydrogen molecule on the Pd(100) 
surface in order to relate it to the reactivity index of Fig. 1. The DFT calculations |28|j 
show that a H2 molecule at a constant height, say at 1.8 A above the center of the top Pd 
layer, experiences a significant corrugation of the potential energy as a function of its lateral 
position (when the H 2 center of mass is at 1.8 A the hydrogen wave functions extend up 
to Z ~ 1.0 A): The potential energy for the H2 center of mass at the hollow site is close to 
zero, at the bridge site it is attractive by 0.12 eV, and at the on top site the attraction is 
0.28 eV (see Ref. [^1). This corrugation can be inferred, at least qualitatively, from the w N 
contour plot of Fig. 1: The strongest decrease of the surface electron density at distances 
between 1.0 and 1.5 A occurs at the on top site. Why does this strongest decrease of w N (r) 
parallel the highest attraction for the molecule ? The dissociative adsorption of H2 involves 
the breaking of the H-H bond and at the same time the formation of new hydrogen-surface 
bonds |T7| , p8| . In a bond-orbital picture jD]] hydrogen dissociation is non-activated if the 



H 2 cr g -orbital as well as the a u interact with the occupied d-states, forming bonding states 
which are occupied (below the d band), and if the corresponding antibonding orbitals remain 
essentially empty, i.e., if their DOSs lie essentially above the Fermi level. In particular the 
partial occupation of the a u state destabilizes the H-H bond. It is quite plausible that this 
scenario is most effective at a site of the metal surface with a high density of d states just 
below the Fermi level, which means at regions with large negative values of w N (r). Indeed, 
the spatial distribution of w N (r) in Fig. |l]d suggests that this region is at the on-top position. 
Here the H 2 molecule interacts strongly with the (3z 2 — r 2 )-orbital of the surface Pd atom 
which gets partially depleted and electrons are transferred into the (x 2 — |/ 2 )-orbitals located 
in the surface plane. 

Figure ^| displays DFT-GGA results for the interacting H 2 -Pd (100) system; shown is the 
electron density change for a geometry with the hydrogen molecule at Z — 1.8 A over the 
on top position. This change is defined as the electron density of the H 2 -Pd(100) system 
minus the density of the clean Pd (100) surface and minus the density of a H 2 molecules 
with the same bond length. It confirms the just-mentioned character of the interaction, in 



particular the decrease of the occupation of the (3z 2 — r 2 )-orbital of the surface Pd atom is 
clearly visible. All features in Fig. |2] are consistent with the spatial variation of w (r) and 
with the interpretation of hydrogen dissociation in terms of a tight-binding picture [|T^] . 

While the reactivity concept is based on perturbation theory and can thus give quan- 
titative information only for weakly interacting systems, it is not unreasonable to expect 
that for strong interactions the concept may still give qualitatively correct guidelines. In- 
deed, we find that for the chemisorption of atomic hydrogen at Pd (100) the function w N can 



explain why the favorable adsorbate site is the fourfold hollow. The DFT calculations |29 
of this chemisorption system tell that the hollow site adsorption is by about 0.3 eV more 
favorable than the bridge and 0.6 eV more favorable than the on-top site adsorption. Fol- 
lowing the same lines as used for the discussion for the H 2 molecule we expect that the H 
atom should prefer a site at which it has good overlap with dark (i.e. negative) sections of 
w N (r). Because of the smallness of the H-atom, this is the hollow site (compare Fig. 1), in 
good agreement with the results of the DFT calculations of the interaction system. Figure |3| 
gives the calculated electron density change induced by a monolayer of H adatoms placed 
in the surface hollow sites. It compares surprisingly well with the prediction of the local 
polarization in Fig. Qa. The comparison shows that the spatial variation of the function 
^^(r) is directly correlated with the shape of the electron polarization introduced by the 
adsorbed hydrogen. 

In conclusion, the function w N (r) [see equ. 2, 3] of the clean surface is introduced 
as an isoelectronic reactivity index supplementing the local softness. It is closely related 
to the local polarizability of valence electrons induced by an atom or molecule chemically 
interacting with the surface. For the example of the dissociative adsorption of H2 molecules 
on Pd (100), we show that the spatial variation of the function w N (r) allows rationalization 
of the preferred dissociation pathways as well as the preferred chemisorption sites of H 
adatoms. The above discussion was done for metal substrates. However, the same reasoning 
applies also for semiconductors and insulators. 
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(a) (b) 




(c) (d) 




FIG. 1. Contour plot of the function w (r) for clean Pd(100) along planes parallel to the 
surface at different heights Z above the center of the surface Pd atoms: a) Z=0 A, b) Z=0.5 A, 
c) Z=l A, and d) Z=1.5 A. Units are lO^A" 3 eV~ 2 . 
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FIG. 2. Electron density change of a H2 molecule above Pd(100) [see text] along a the (Oil) 
plane. The molecule is placed at a height Z = 1.8 A above a surface Pd atom with the molecular 
axis parallel to the surface, and the two H atoms point toward hollow sites. Units are 10 _3 /A 3 . 
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FIG. 3. Electron-density change at Pd (100) due to the adsorption of one monolayer of hydrogen 
at the surface hollow sites. The plane is the same as in Fig. la, the units are 10 _3 /A 3 . 
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